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^Although  methods  of  filtering  have  been  developed  for  representing  on  a 
planar  surface  the  sub-grid  scale  process  of  diffusion  in  numerical  modeling 
of  the  atmosphere,  the  proper  form  of  filtering  for  a  spherical  domain  remains 
to  be  selected.  In  this  study,  the  characteristics  of  several  different  forms  of 
a  highly  scale  dependent  low  pass  filter  are  examined  and  compared  for  the 
case  where  the  filter  is  applied  to  a  scalar  field  on  the  surface  of  a  sphere. 

The  phase  and  amplitude  response  functions  of  the  various  forms  of  the  filter  - 
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Abstract  (Continued) 


icate  that  the  simplest  form,  although  it  does  not  preserve  area -weighted 
mean  values,  approaches  most  closely  the  criteria  established  for  the  Ideal 
filter.  This  indication  is  verified  by  test  computations  in  which  each  form 
of  the  filter  Is  applied  up  to  10,  000  times  to  a  noisy  scalar  field. 
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Linear  Filtering  on  the  Surface 

of  a  Sphere 


1.  INTRODUCTION 

An  efnc  lent  and  uaeful  niter  haa  been  developed  to  represent  the  sub-grid 
scale  process  of  diffusion  In  numerical  models  of  the  atmosphere  (Shapiro1).  The 
filter  has  the  property  of  variable-scale  dependency.  That  la,  as  the  order  of  the 
operator  Is  Increased,  the  dissipation  of  the  short  wavelength  components  la  In¬ 
creased  relative  to  that  of  the  longer  wavelengths.  The  operator  of  order  8  effec¬ 
tively  eliminates  components  with  wavelengths  less  than  4-grld  Intervals  but  pro¬ 
duces  no  measurable  damping  of  components  with  wavelengths  greater  than  6-grtd 
Intervals  (Shapiro3).  When  used  In  numerical  models  of  the  atmosphere,  the 
filter  promotes  computational  stability  by  removing  the  highly  noise  contaminated 
short  waves  (particularly  the  2 -grid -Interval  wave)  but  because  of  Its  scale  depen¬ 
dency  does  not  damp  the  longer,  meteorologically  Important  wave  components 
(Hunt,  3  Francis,  *  Voice  and  Hunt,  ®  Tapp  and  White,  6  Mudrlck,  *  Kalnay-Rtvas 
et  al9).  Consequently,  use  of  the  filter  minimizes  the  unfortunate  tendency  of  all 
stable  models  to  suppress  the  growth  of  barocllnlc  Instability  and  eddy-klnetlc 
energy. 


(Received  for  publication  29  October  1979) 

Due  to  the  number  Qf  references  to  be  Included  as  footnotes  on  this  page,  the  reader 
Is  referred  to  the  list  of  references,  page  33. 
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There  is  no  ambiguity  in  the  use  of  the  Alter  in  a  grid  system  which  is  mapped 
on  a  tangent  plane  or  on  polar  stereographic  coordinates.  The  proper  form  of  the 
filter  is  developed  from  the  linear  finite  difference  analogue  of  the  second  deriva¬ 
tive  in  space  or  time.  However  it  is  not  at  all  clear  what  the  proper  form  of  the 
filter  is  in  spherical  geometry,  because  of  both  the  convergence  of  the  meridians 
and  the  potential  mathematical  singularity  at  the  poles.  The  purpose  of  this  study 
is  to  examine  a  variety  of  possible  filter  forms  in  order  to  derive  a  rational  basis 
for  selection  of  a  best  form. 

The  one -dimensional  operator  of  order  p  =  1,  2,  3,  ...  given  by 


(1) 


o 

where  6  fj  is  represented  as  a  finite  difference  analogue  of  the  second  derivative 

of  f  (in  space  or  time)  at  the  gridpolnt  i  multiplied  by  the  square  of  the  space  or 

time  increment.  The  above  operator  has  been  shown  to  be  a  useful  scale -dependent 
o 

filter  (Shapiro  ).  For  any  integer  p>  1  the  filter  removes  2  grid -interval  compo¬ 
nents  completely  and  becomes  increasingly  scale  dependent  as  p  increases.  With 

o 

a  suitable  choice  of  the  form  of  6  ,  the  Alter  preserves  the  phase  of  all  Fourier 
components,  introduces  no  new  components  nor  does  it  amplify  any  existing  com¬ 
ponent.  It  is  therefore  highly  effective  as  a  dissipative  and  diffusive  mechanism 
(Shapiro1).  It  is  particularly  useful  in  any  highly  repetitive  or  iterated  procedure 

such  as  the  numerical  marching  procedures  of  numerical  weather  prediction  (NWP) 

6 

models.  Although  there  is  no  confusion  concerning  the  appropriate  form  of  6  in 

one  dimension  with  uniform  gridspacing  or  in  more  than  one  dimension  in  Cartesian 

2 

coordinates,  there  are  a  variety  of  forma  that  6  may  assume  when  applied  on  the 
surface  of  a  sphere. 

4  2 

Francis  has  applied  the  filter  in  a  form  in  which  6  is  an  analogue  of  the 

finite  difference  representation  of  the  second  derivative  in  the  north-south  direction 
on  a  sphere.  This  form  has  the  property  of  conserving  area-weighted  mean  quan¬ 
tities  and  with  a  suitable  choice  of  grid  location  and  spacing  it  is  self-consistently 

o 

well-behaved  at  the  poles.  Kalnay-Rtvas,  et  al  have  compared  the  above  form  of 
the  filter  with  the  Cartesian  form  which  neglects  the  spherical  geometry  and  a 
simpler,  but  not  consistent  version  of  the  Bpherical  form.  They  found  little  dlf- 

A 

ference  among  the  three  forms  although  only  the  form  suggested  by  Francis  is 
consistent  with  the  spherical  geometry. 

2 

Since  the  properties  of  the  filter  are  largely  determined  by  the  form  of  6 

2 

this  paper  explores  the  characteristics  of  the  filter  with  various  forms  of  6  as 
applied  to  operations  on  the  surface  of  a  sphere. 
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2.  FORMS  OF  6  2 


If  the  function  g(x,  y,  z)  in  rectangular  Cartesian  coordinates  is  expressed  in 
spherical  coordinates  in  terms  of  t(9,  X,  r)  we  have 


x  =  r  cos  $  cos  X 
y  =  r  cos  9  sin  X 
z  =  r  sin  $ 


(2) 


where  r  is  the  distance  from  the  origin  normal  to  the  surface  of  the  sphere,  f)  is 
the  latitude  measured  positive  northward  from  the  equator  and  X  is  the  longitude, 
measured  positive  toward  the  east.  The  LaPlacian  of  g  or  f  is  given  by 


v2g 


O.-r. 

\3x2  9y  02.  / 


where 


v2f 


.  <r2  c2  $  *  <-2  co.  .)-*  Is  (-  .  |l)  ♦  r-2  I;  (c2  ft) 


<3> 


Let  fy  represent  the  function  f  on  the  surface  of  a  sphere  at  the  gridpolnt  (ij) 

where  i  is  an  index  for  X  and  j  an  index  for  9.  AX  and  AO,  the  grid  Increments  in 

the  i  and  j  directions,  are  constants  in  space  though  not  necessarily  equal  to  each 

other.  The  first  and  second  terms  on  the  right  of  (3)  represent  respectively  the 

variation  of  f  in  the  east-west  direction  along  a  latitude  circle  and  in  the  north- 

south  direction  along  a  meridian.  If  we  restrict  attention  to  the  surface  of  a  sphere, 

the  third  term  is  zero.  In  applying  the  filter  (1)  to  the  field  of  f..  in  the  east-west 

2 

direction,  there  is  no  ambiguity  concerning  the  choice  of  6  .  We  have 


flj  =  fi+l.j  '  2fij  +  fi-l.j 


(4) 


where  the  subscript  to  the  left  of  6 2  indicates  the  direction  of  the  filtering.  We 

2 

have  examined  the  following  forme  of  f^ 


Form  1 


fi,j+l  ‘  2ftj  +  ft.j-l 
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I 

i 


Form  2 


[(fi,  j+1  '  V  cos  Vl/2  '  (flj  "  fi,  j-l*  COS 


Form  3 

zzhr.  [fi.  j+i  008  V»  -  2fij 008  6j +  fi.  j-i coa  V] 

Form  4 


i.  j+1 


-2fiJ  +  fU-i 


A0 


T<fi.j-n 


tan  e. 


Form  5 


fi,  j+1  -  2fij  +  fi.  j-1  +  (fi,  j+1  -  fi.  j-l)(cos  6 j+1  '  cos  Vl)/4  cos  flj 

Form  1  is  the  Cartesian  form  and  avoids  singularities  at  the  poles  by  neglecting 
the  spherical  geometry.  Forms  2.  4  and  5  are  finite  difference  versions  of  the 
second  term  on  the  right  of  (3).  Form  3  is  an  arbitrary  simplification  of  Form  2. 
Forms  2  to  5  all  have  potential  singularities  at  the  poles,  but  if  the  gridpoints  j 
are  chosen  so  that  the  north  and  south  poles  are  situated  at  1/2-grid  interval  north 
and  south  respectively  of  the  most  northerly  or  southerly  gridpoint  j,  then  Form  2 
used  by  Francis,  4  has  no  singularity  and  may  be  applied  in  a  self-consistent  fashion. 
Some  arbitrariness  is  unavoidable  near  the  poles  with  the  use  of  the  other  forms. 


3.  PHASE  AND  AMPLITUDE  RESPONSE 

o 

The  order  p  amplitude  response  of  the  filter  (1)  with  the  use  of  Form  1  for  6 
has  been  shown  (Shapiro1)  to  be 

Rj  =  1  -  sin2po  (5) 

where  a  =  nA/2,  n  is  the  wave  number  of  the  Fourier  component  whose  wavelength 
is  given  by  X  -  2*/n  and  A  is  the  gridspacing.  Equation  (5)  shows  that  no  Fourier 
component  can  be  amplified  by  the  filter,  that  the  wave  of  length  2A  is  removed 
regardless  of  the  value  of  p  and  that  as  p  increases  the  highest  wave  number  for 
which  there  is  virtually  no  damping  also  increases.  Therefore  the  filter  becomes 
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increasingly  scale  dependent  as  p  increases.  That  is,  the  maximum  slope  of  the 
amplitude  response  function  becomes  increasingly  steep,  thus  effectively  distin¬ 
guishing  between  those  frequencies  which  suffer  virtually  no  damping  from  those 
which  are  almost  completely  damped.  The  phase  response  or  phase  shift  of  this 
form  of  the  filter  is  zero  since  this  form  is  symmetrical.  This  Important  property 
ensures  that  no  new  Fourier  components  are  Introduced  merely  by  filtering.  These 
desirable  properties  of  the  filter  do  not  obtain  with  any  of  the  other  forms  of  4^. 

The  phase  ($)  and  amplitude  response  (R)  functions  for  the  filter  (1)  with  Forms 
2 

2  to  5  of  6  and  with  p  =  1  are  given  below 


Form  2s 


=  tan 


-1 


sine  cos  a  sin  ten  0^ 

wr 


1  -  sin  a  cos 


(6) 


r 

2  f 

1  -  sin  a  cos  1 

•t) 

^  2  2  o  // 

+  sin  a  cos  a  sin  { = 

rj  tan2  J 

\ 

V 

1/2 


Form  3: 


=  tan" 


sin  2a  sin  A0  tan  0, 


i  +  cos  2a  cos  A0 


Rg  |^[1  +  cos  2a  cos  A0l  ^  +  sin**  2a  sln^  A0  tan**  0.^ 


1/2 


Form  4: 


*4  =  tan'1  ^^2  tan  a  tan  0^ 


4  _  A  _i_2 _ 2 _ 2 


R4  “  cos  a  +  sin  a  cos  a  tan  0. 


9j  M*] 


1/2 


(7) 


(8) 


(9) 


(10) 


(11) 


11 


Form  5: 


4.  *  tan 


1  ^  j  sin  40  tan  a  tan 

(4  1  2  2  2  2  1 

cos’  a  +  j  sir  a  cos  a  tan  0^  sin  40  J 


(12) 


(13) 


Although  in  the  limit  as  a  approaches  zero  (and  in  the  case  of  Form  3  as  40 
also  approaches  zero)  the  forms  of  the  filter  (Forms  2  to  5)  have  the  same  ampli¬ 
tude  response  and  phase  shift  as  Form  1,  in  each  case  for  non-zero  values  of  a  and 
40  the  phase  shift  is  different  from  zero  everywhere  except  at  0^  =  0.  Further¬ 
more,  the  phase  shift  is  a  function  of  latitude  and  wave  number.  This  means  that 
application  of  the  filter  (Forms  2  to  5)  to  a  function  with  a  single  Fourier  compo¬ 
nent  will  introduce  new  Fourier  components  merely  by  the  filtering  process.  In 
addition  the  amplitude  response  is  not  only  a  function  of  wave  number,  as  it  should 
be,  but  also  a  function  of  latitude.  Beyond  this,  and  potentially  more  serious  is 
the  fact  that  the  amplitude  response  of  the  filter  (with  Forms  2  to  5)  is  a  function 
of  tan  0.  which  becomes  increasingly  large  as  0^  approaches  ±ir/2.  This  feature 
of  Forms  2  to  5  introduces  the  possibility  that  the  filter  could  actually  amplify 
certain  Fourier  components. 

If  we  assume  that  0,  is  chosen  so  that  it  does  not  include  the  poles  (0  =  ±ir/2), 

J  o 

then  the  maximum  value  of  tan  0.  (and  therefore  tne  maximum  value  of  R  for 

Forms  2  to  5)  occurs  when  0.  =  ±1/2  <*  -  40).  That  is,  for  these  Forms  R  is  a 

3  2  2 

maximum  for  0^  1/2-grid  interval  from  the  poles.  But  tan  [±l/2  (jr  -  40))  =  ctn 

(40/2).  Therefore,  the  maximum  value  of  (7) 


2  2 

Since  0  £  cos  (40/2)  (2  -  cos  (40/2)1  £  1,  Rglmax)  £  1  and  approaches  Rj  as 
(40/2)  -*0.  Nevertheless,  Rg  is  a  function  of  0^  and  this  property  as  well  as  the 
phase  shift  that  occurs  with  the  Form  2  filter  will  introduce  false,  extraneous 
Fourier  components.  The  fact  that  Rg  is  well-behaved  in  terms  of  not  amplifying 
any  existing  component  implies  that  the  extraneous  components  will  remain  small. 
This  feature  does  not  obtain  with  Rg  (max)  which  is  given  by 


Rg  (max)  = 


a  +  cos 


2 


(15) 


12 


For  reasonable  values  of  AS,  (for  example,  AS  <  11°),  cos2  (AS/2)  >  0. 99.  With 
cos2  < AS/2)  =  0. 99,  Rg  (max)  >  1  for  0. 0102  <  sin2  a  <  0. 6587,  which  Includes  the 
bulk  of  the  Fourier  components  of  Interest. 

Rj  (max)  Is  given  by 


(16) 


Inasmuch  as  the  quantity  0  <  (AS/2)2  cos2  (AS/2)/sin2  (AS/2)  si,  it  is  apparent 
2 

that  R^  (max)  s  1.  Similarly,  the  closely  related  Rg  (max)  which  is  given  by 


R&  (max)  = 


(1  -  sin  a)  1  -  sin 


(17) 


2 

can  easily  be  seen  to  satisfy  the  condition  R&  (max)  s  1. 

It  appears  that  although  all  forms  of  the  filter  with  the  possible  exception  of 
Form  3  can  be  expected  to  behave  reasonably  well,  only  Form  1  exhibits  no  poten¬ 
tially  pathological  behavior.  It  is  of  interest  therefore  to  compare  the  various 
forms  of  the  filter  in  severe,  but  realistic  circumstances. 


4.  TEST  COMPARISONS 

For  purposes  of  test  comparison  we  make  use  of  a  convenient  scalar  function 
of  latitude  (ft)  and  longitude  (X)  which  had  been  used  by  Yee10  in  another  study. 
This  function  has  the  form 


f(X,  0)  =  sin  0 


18 

7:  B(m)  cosm 
m  =  l 


0  cos  ml  [1  +  o(m)  +  3(0)1 


(18) 


where  A  is  a  constant,  B  is  a  weighting  factor  which  is  a  function  of  m,  the  wave 
number  parameter,  and  a  and  3  are  random  phase  shifts.  f(X,  0)  was  evaluated 
with  a  resolution  of  A0  -  AX  =  2ir/36,  yielding  values  for  f„  for  i  *  1,  2, . ...  36 
corresponding  to  X  =  0,  10 . 350  deg  and  j  =  1,  2 . 18,  19 . 36  correspond¬ 

ing  to  0  =  -85,  -75, ....  85,  85, ....  -85  deg.  Thus  gridpoints  (i,  j),  for  which 


10.  Yee,  S.  Y.  K.  (1978)  An  Efficient  Barotropic  Vorttcity  Equation  Model  on  a 
Sphere.  Env.  Res.  Pap.  No.  645,  AFuL-f R-78-0273,  AD  A065  117. 
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1  s  j  «  18,  are  on  the  1th  meridian,  but  grldpoints  for  which  19  s  j  js  36  are  on 
the  (i  +  18)th  meridian.  Together,  theae  grldpotnta  form  a  complete  meridian 
circle,  fy  is  therefore  doubly  periodic  with  l  =  t  ±  36  and  j  =  j  ±  36.  Since  there 
are  only  648  (18  times  36)  distinct  values  of  fy  each  distinct  value  is  assigned 
twice.  That  is  fy  =  fl+lg  37_j.  fy  was  treated  in  this  fashion  so  as  to  permit  the 
calculation  of  the  double  Fourier  spectrum.  The  field  fy  (j  *  1,  2, . . .,  18)  is  shown 
in  Figure  1.  fy  has  been  scaled  so  that  -1  <  fy  <  1,  although  for  our  purposes 
the  numerical  value  of  fy  is  less  important  than  its  spatial  or  spectral  distribution. 
It  can  be  seen  from  Figure  1  that  f  is  "noisy"  with  abrupt  variations  between  rela¬ 
tive  maxima  and  minima.  This  "noisy"  appearance  is  consistent  with  its  two- 
dimensional  Fourier  spectrum  shown  in  Figure  2  which  indicates  that  an  apprec¬ 
iable  fraction  of  the  total  variance  of  fy  is  contained  in  wave  numbers  9  to  18 
corresponding  to  wavelengths  from  4-  to  2-grid  intervals. 

We  also  make  use  of  a  modified  function,  f*(V,  6),  given  by 

f*(X..  0)  =  t(\,  0)  cos  0  (19) 


which  is  illustrated  in  Figure  3  and  whose  two-dimensional  Fourier  spectrum  is 

shown  in  Figure  4.  The  principal  difference  between  ff.  and  f..  is  the  much  smaller 
*  lJ  lJ  * 

variation  of  fy  with  respect  to  longitude  near  polar  latitudes.  In  this  regard  fy  is 

more  consistent  with  the  spherical  geometry.  Because  of  the  relatively  large 

$ 

amount  of  variance  at  high  wave  numbers,  both  fy  and  fy  represent  severe  tests 
of  the  effectiveness  of  the  filter. 

*  2 

In  applying  the  filter  to  fy  and  fy  we  use  the  form  of  6  given  by  (4)  in  the 
i-direction  regardless  of  the  form  used  in  the  j-direction.  In  general,  the  one¬ 
dimensional  operator  of  order  p,  Eq.  (1),  can  be  applied  as  a  sequence  of  p, 

3-point  operations  or  as  a  single  (2p  +  1)  -  point  operation.  Because  of  the  special 
treatment  required  near  the  poles  it  was  convenient  to  apply  the  operation  in  both 

the  i  and  j  directions  as  sequences  of  3 -point  operations.  In  order  to  evaluate 
o 

.6  f.,,  it  is  necessary  to  have  available  values  at  f.  .  .  and  f.  as  well  as  at  f... 

J  *  J  J  +  1  l»  J"4 

When  j  =  1  or  18,  f.  .  ,  and  f.  .  .  respectively  are  considered  to  be  the  gridpoints 

across  the  poles  situated  at  f.+jg  j.  This  expedient,  which  is  not  necessary  in  the 

case  of  Form  2,  was  adopted  primarily  because  it  seems  reasonable  to  treat  the 

poles,  insofar  as  possible,  in  the  same  fashion  as  any  other  point.  However,  some 

tests  were  carried  out  in  which  the  poles  were  treated  as  quasi-solid  boundaries 

with  generally  unsatisfactory  results. 

Two-dimensional  filtering  may  consist  of  different  orders  of  filtering  in  each 
direction,  consistent  with  the  two-dimensional  Fourier  spectrum  of  the  field  being 
operated  upon.  However,  for  the  present  study  the  same  order  of  filtering  (with 
p  =  8)  was  used  for  all  fy  for  both  the  i  and  j  directions.  Filtering  was  first 
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Figure  2,  Two-Dimensional  Fourier  Spectrum  of  Initial  Field,  fy, 
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Figure  3.  Unsmoothed,  Modified  Stream  Function  Field,  f. 
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Figure  4.  Two-Dimensional  Fourier  Spectrum  of  Modified  Field 
r..  Multiplied  by  100 


completed  In  the  i-dlrectlon  and  then  these  partly  filtered  values  were  filtered  in 
the  j  -direct  using  one  of  the  5  Forms  listed  in  Section  2. 

For  simplicity  the  two-dimensional  filtering  operation  making  use  of  Eq.  (4) 

2  2 

for  i6  and  Form  (  )  for  j6  will  be  referred  to  as  a  Form  (  )  filter. 

Figures  5  to  7  and  8  to  10  show  the  result  of  applying  the  Form  1,  2  and  3 
filters  (to  fy)  1000  and  10, 000  times  respectively.  Examination  of  Figures  5  to  7 
shows  a  remarkable  change  brought  about  by  the  filtering  as  compared  with  the 
initial  field  (Figure  1).  The  filtered  fields  are  smooth.  Furthermore,  the  differ¬ 
ences  among  the  three  forms  are  small,  except  for  Form  3  at  high  latitudes.  The 
misbehavior  of  Form  3  at  high  latitudes  was  anticipated  by  Eq.  (15).  Comparison 
of  Figures  5  to  7  with  Figures  8  to  10  shows  that  after  1000  applications  of  any  of 
the  three  forms  little  further  filtering  takes  place. 

Figures  11  to  16  show  the  two-dimensional  Fourier  spectra  of  the  filtered 
results  shown  in  Figures  5  to  10  respectively.  Figures  11  to  16  confirm  the  con¬ 
clusions  derived  from  Figures  5  to  10,  but  in  addition  they  clearly  reveal  other 
characteristics  of  various  forms  of  the  filter. 

For  our  purposes,  an  ideal  filter  would  remove  high  wave  number  components 
both  in  the  i-  and  j-dlrections,  would  preserve  the  phase  and  amplitude  of  lower 
wave  number  components  and  would  not  create  any  components  which  were  not 
present  before  filtering.  Although  all  three  forms  of  the  filter  are  highly  effective 
in  removing  the  high  wave  number  components  in  the  i  (east -west)  direction,  only 
Form  1  is  equally  effective  in  the  j -direction.  Furthermore  only  Form  1  succeeds 
in  avoiding  the  amplification  of  existing  components  or  introducing  new  components, 
although  Form  2  is  obviously  superior  to  Form  3. 

Figures  17  to  22  show  the  result  of  filtering  the  f!*  field  (Figure  3)  with  the 
Form  1,  2,  and  3  filters.  These  figures  are  comparable  to  Figures  5  to  10 
respectively  excqst  for  the  cos  0  modification  of  the  initial  field.  The  comments 
made  concerning  Figures  5  to  10  apply  with  equal  force  to  Figures  17  to  22.  Judg¬ 
ing  from  these  figures,  there  would  appear  to  be  little  basis  for  choosing  among 
the  3  forms  of  the  filter  except  for  the  simplicity  of  Form  1.  However,  Figures 
23  to  25  and  Figures  26  to  28  which  show  the  two -dimens  tonal  Fourier  spectra  of 
the  fjj  fields  expressed  as  percentages  of  the  ortglnal  spectrum  after  1000  and 
10,  000  applications  of  the  Form  1,  2  and  3  filters  demonstrate  the  superiority  of 
the  Form  1  filter.  This  form  of  the  filter  appears  to  have  the  characteristics  of 
the  ideal  filter  outlined  above. 

Figures  29  to  32  show  the  comparable  two-dimensional  Fourier  spectra  of  the 
fjj  fields  after  filtering  with  Forms  4  and  5  respectively.  Forms  4  and  5  appear  to 
have  characteristics  which  are  similar  to  those  of  Form  2.  This  is  not  surprising 
in  view  of  the  similarity  of  their  phase  and  amplitude  response  functions. 
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Figure  6.  f{i  After  1000  Applications  of  the  Form  2  Filter 
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Figure  11.  Two-Dimensional  Fourier  Spectrum  (X  100)  of  f..  After 
1000  Applications  of  the  Form  1  Filter  y 
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Figure  12.  Two-Dimensional  Fourier  Spectrum  (X  100)  of  f..  After 
1000  Applications  of  the  Form  2  Filter 
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Figure  13.  Two -Dimena lonal  Fourier  Spectrum  (X  100)  of  f,,  After 
1000  Applications  of  the  Form  3  Filter 
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Figure  29.  Percent  of  Initial  Fourier  Amplitudes,  Remaining 
After  1000  Applications  of  Form  4  Filter  J 
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Figure  30.  Percent  of  Initial  Fourier  Amplitudes,  f!*..  Remaining 
After  1000  Applications  of  Form  5  Filter 
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Figure  31.  Percent  of  Initial  Fourier  Amplitudes,  f. .,  remaining 
After  10,  000  Smoothings  of  Form  4  Filter 
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5.  CONCLUSIONS 


Among  the  five  forms  of  the  filter  which  were  tested  on  noisy  fields  on  the 
surface  of  a  sphere.  Form  1,  the  simplest  form,  stands  out  as  superior.  Form  S, 
which  has  a  tendency  to  pathological  behavior  near  the  poles  is  inferior  and  should 
be  avoided.  Forms  2,  4  and  5  appear  to  be  quite  similar  in  their  effects. 

In  applications,  the  order  of  the  filler,  p,  should  be  chosen  so  that  there  is  a 
sharp  transition  between  those  Fourier  components  which  are  considered  undesir¬ 
able  (and  ace  strongly  damped)  and  those  components  which  one  wishes  to  remain 
unaffected  by  the  filtering  process. 
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